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ABSTRACT 


Two-dimensional turbulent wakes in constant adverse and 
favorable pressure gradients were studied experimentally. 

The static and dynamic pressures through the wake were ob- 
tained with calibrated Prandtl-type pitot-static probes, 

and a computer program was written to extract from these 
data the velocity profiles, displacement thickness, momentum 
thickness, and the shape factor. This program further pro- 
vides the turbulent shear coefficient, eddy viscosity, and 
Prandtl mixing length within the wake. 

A wind tunnel was designed specifically for this experi- 
ment and its construction details are included. The tunnel 
Beet section is unique in that it can be adjusted to generate 
any specified adverse or favorable pressure gradient. 

The results of the experiment indicated that the eddy 
viscosity was a function of the history of the flow rather 
than of the local eon flow field, as is clessicalayeas SÚmeo . 
Bie reaction of the turbulent wake to a favorable pressure 
gradient Are encountering an adverse pressure gradient was 
Significantly different from the reaction of the wake to 


only a favorable pressure gradient. 
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NOMENCLATURE 


cross-sectional area 
indicates absolute value 
turbulent exchange coefficient (Eq. 3, Introduction) 


constant in Prandtl's eddy viscosity (Eq. 5, Intro- 
Guctiion ) 


width of mixing region (Eq 5, Intredueticn) 


constant in test section area distribution 
(Eq. 7, Appendix A) 


volumetric flow rate through the test section 


negative of the test section pressure gradient 
(Appendix A) 
turbulent shear coefficient = r/5pu“e 


wake generator diameter 
indicates differentiation with respect to x 


indicates partial differentiation with respect to x 
displacement thickness 

eddy viscosity = 1/38 

indicates exponentiation 

arbitrary example function 

shape factor = 6%*/6 

indicates i direction (i = 1,2,3) 

large distance in the cross-streamwise direction 
indWestes j' déeection (j 2 1,2,3) 

length of pressure gradient 


Prandtl's mixing length (Eq. 4, Introduction) 
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logarithm to base e 
laminar eocetftfielent of viscosity 
wake centerline conditions 
static pressure 
à Еа 
dynamic pressure = 50V 
fluid density 
stress tensor in Navier-Stokes equations 


mean flow stress tensor in the time averaged 
Navier-Stokes equations 


indicates time 


indreateserkescharacteristic period used for 
time averaging 


Reynolds shear stress tensor 

momentum thickness 

mean velocity in the x-diréerion 

fluctuating velocity in the x-direction 

mean velocity in i,j directions (i,j = 1,2,3) 
fluctuating velocity in 1, directions (1. το. 
mean velocity in the y-direction 

fluctuating velocity in the y-direction 

one dimensional mean flow velocity 

streamwise coordinate 

j coordinate 

cross-streamwise coordinate 


lateral position at which C. is being calculated 
(Eq. 5, Appendix B) 





Subscripts 


e = local external flow conditions 

i,j = directionality parameter (i,j = 1,2,3) 

max = maximum value within the wake 

min = minimum value within the wake 

T Menear estate cq aras, 

x,y = indicates partial differentiation with respect to 
x and y 

1 = test section inlet conditions 

2 = test section exit conditions 

superscripts 


indicates a time averaged quantity 


designates a fluctuating quantity 





I. INTRODUCTION 


When tangential stresses are applied to a fluid having 
@ivernal friction, shearing motions are set up in line with 
спе stresses and in conformity to the shape of the bound- 
aries. Within this flow field, various kinds of secondary 
mocions become possible. Regular ones, when they occur, 
are readily accounted for in terms of pressure gradients 
resulting from the curvature of the main flow, usually 
imposed by the shape of the boundaries. Irregular ones, 
med turbulence, are by far the more common, and their 
ject cause is less obvious. Their occurrence does nov 
depend on the shpae of the boundaries, but like all secon- 
dary flows they must depend on a generating mechanism within 
the flow itself. According to present evidence the initial 
onset of turbulence occurs suddenly by a breakdown of the 
laminar flow in localized regions. Turbulence is convected 
downstream in the manner of any other fluid property and 
washed away from its point of origin. Repeated breakdowns 
are ШОО Сох general, required to maintain a continuous 
supply of turbulence, and instability of laminar flow is 
em essential part of this process. 

In treating turbulence mathematically it is generally 
assumed that the motion can be separated into a mean flow, 


whose components are u, (Cartesian tensor notation), and a 


3t 


superimposed fluctuating motion whose components are uj, the 


al 





mean values of which are zero. For the flow of an incom- 
peessible Newtonian fluid, this. decomposition»can be applied 
to the dependent variables of the Navier-Stokes equations 
and time-averaged, the latter defined as 

T 1 CFT 

(ы m f Horti CIL) 

t 

where T is a period of time large compared with the char- 
meveristic time of the turbulent oscillations. “fhe resulting 


time-smoothed Navier-Stokes equations are 
p = 2 (o.. - pu.u, - pulu! (2) 
ot ðX ag de] d 


These equations have the same form as the original 





Navier-Stokes equations if ZE is replaced by 914 - pu,u,. 


Thus the equations for the mean flow are the same as the 


original equations of motion, except that there are addi- 


A e саесатр 


tional virtual stresses, T = ur which represent the 


1] 
mean rate of transfer of momentum across a surface due to 
[ee velocity fluctuations. These virtual stresses were 
first introduced by Osborne Reynolds and are known by his 
name [Ref. 4]. 

i 


in the i- and j- directions at one point in space and uu; 


The quantities u! and u; are the velocity fluctuations 


is a measure of the correlation between the fluctuations. 
Correlations between fluctuations at two points in space 


are also measurable and give the size and orientation of 
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me turbulent eddies. , In trying to solve analytical ter 
Exe Second-order correlation tensor, the solution is given 
in terms of third-order correlations. Simiderivyvsguhe 
Solution for the third-order correlation is given in terms 
of the fourth-order correlation, and so on. Hence the 
system cannot be closed and it is impossible to solve 
exactly for the second-order correlations [Ref. 1]. 

The various methods for the*eadleulationwef спе ent 
flow field developed to date are based on empirical hypoth- 
eses. These approaches endeavor to establish a relationship 
between the Reynolds stresses produced by the fluctuating 
velocities and the mean values of the velocity. components 
end their spatial derivatives. Ву this means, the virtual 
meresses are given an analytical form which closes the 
system of equations and constitutes the starting point for 
Bue calculation of the mean flow. 

The first hypothesis relating the Reynolds stress to 
mean flow properties was advanced by Boussinesq and intro- 
sed a turbulent mixing coefficient, А, such that the 


turbulent shear stress could be written 


i du 
q = A, dy (3) 


This coefficient corresponds to the viscosity, u, in laminar 
flow and is therefore often called the "virtual" viscosity. 


Equation (3) has the great disadvantage that the virtual 
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pscosity, A. is not a property of the fluid, аве ADe 
depends on the structure of the flow field. 

In 1925 Prandtl proposed his mixing length hypothesis 
in which the turbulent stresses are modelled as follows: 

t= p° DIE (4) 
The mixing length, X, is that distance in the transverse 
direction which must be covered by an aggregation of fluid 
mercicles, travelling with its original mean velocity, in 
order to make the difference between its velocity and the 
velocity in the new lamina equal to the mean transverse 
muctcuations in the turbulent flow. The mixing length is 
meme not a property of the fluid ΡΞ ΙΙ; ΠΠ οπου, 
eppurely local function and it is far simpler to make 
plausible assumptions concerning the mixing length, X, than 
mor the virtual viscosity, Ат. This constitutes the essen- 
tial superiority of Prandtl's equation over that due to 
Boussinesq. Even so, Equation (4) is unsatisfactory in that 
the apparent kinematic viscosity, Е = gu vanishes at 
d 
d 


maximum and minimum velocity, which is contradictory to 


points where Fra is equal to zero. that 16 20 polnve "oe 


experimental evidence. A simpler expression for e, also 


Proposed by Prandtl which overcomes this problem is 


€ bU. - T) (5) 


1:2 





and is valid only for free turbulence, where a is an experi- 
Memcalily determined, dimensionless constant and b is the width 
Buche mixing zone. Thus €, for all the classical models or 
Reynolds stresses, remains constant over the whole width of 
meee wake and is a function of the local mean flow. 

Several other similar hypotheses have been proposed as 
analytical expressions for the Reynolds stresses which will 
only be mentioned here. They include: Taylor's vorticity 
Meensport theory, Von Karman's similarity hypothesis, and 
Pees universal velocity distribution law, each of which is 
inked to the mixing length concept due to Prandtl. Each 
of the aforementioned theories is valid only for free turbu- 
lut flow. While none of the theories adequately models the 
ОНО ретп process in general, they are useful tools in the 
integration of the equations of motion in special cases. 

In this investigation it was desired to ascertain, at 
Beast qualitatively, whether the virtual viscosity could be 
Gemsidered to be a function of the local mean flow propert es 
in the case of two-dimensional turbulent wakes in constant 
favorable and adverse pressure gradients. A wind tunnel 
was designed and fabricated for this study which incorporated 
a variable-geometry test section that could be set to provide 
any specified pressure gradient. While this study was carried 
out using pitot-static probes as data sensors, it is antici- 
pated that the flexibility and utility of the wind tunnel 
test section will allow a very detailed analysis of the 


turbulent wake using a hot-wire anemometer. This will also 


їй 





permit the assessment of more modern eddy viscosity models 
which involve hypotheses with regard to the dependence of 


the shear on the turbulent kinetic energy. 
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II. THE WIND TUNNEL 


A. INLET 

iieeanlet betimouth shown in Fie. 1 wae od Cpu ECC 
from eight-inch diameter aluminum pipe halves formed into a 
two-feet square inlet section. The corner seams, which were 
K C Sed to the inlet air flow, were filled with epoxy putty, 
hand formed, and sanded into a smooth fillet. This was done 
to smoothly transition the flow around these corners into 
parallel flow. 

Immediately downstream of the inlet bellmouth was the 
plenum section of the tunnel. It consisted of four two-feet 
Square by one foot long steel sections bolted together. 

Cork gaskets were inserted between the sections to prevent 
the influx of air at the joints and allow for smooth air 
flow over these junctions. Finally, two copper screens of 
#16 mesh between the first two and last two sections m enueia 
the ambient turbulence level in the flow prior to its 
entering the nozzle. 

Initial operation of the tunnel indicated the presence 
swirling in the air flow through the test section. This 
swirling was eliminated by the insertion of a four-inch 
thick, two-feet square section of quarter-inch honeycomb 


met prior to the first copper screen. 
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B. NOZZLES 

Two different nozzles were constructed and used during 
the course of this investigation. One nozzle, used with the 
To rable pressure gradient test section, was constructed 
En half-inch Plexiglas.: The nozzle, shown in Pig. 29s 
two feet long, has an exit area of 168 square inches, (24 
inches high by 7 inches wide), and a contraction ratio of 
3.42. The top and bottom sides of the nozzle are parallel 
While the vertical sides contract to form the nozzle. The 
wake generator was attached to the nozzle in its exit plane 
e” axial centerline of the vertical walls. static 
Pressure ports were located at the nozzle entrance and exit 
pues along the axial centerline of the bottom wall. These 
were used to determine the nozzle exit (test section entrance) 
dynamic pressure, which determines the magnitude of the test 
section pressure gradient (see Appendix A). 

Ehe second nozzle was made Trom steel ana Wao Usea wits 
an adverse pressure gradient test section configuration. 
MN nozzle, shown in Fig. 1, is two feet long, has an exit 
area of forty nine square inches (seven inches square) and 
a contraction ratio of 11.76. Static pressure ports were 
Meeated on the right side of the nozzle at the inlet and 
exit for the purpose of determining the exit dynamic pressure. 
The wake generator was not mounted in the nozzle because of 
the difficulty encountered in trying to center and align it 
in the metal; rather in this case it was mounted just at 


the entrance of the test section. 
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C. WAKE GENERATORS 

The wake generators were smooth Plexiglas cylinders 
spanning the width of the tunnel. A one inch diameter 
cylinder was used with the steel nozzle and a two inch 
diameter cylinder was used with the Plexiglas nozzle; thus 
"equivalent" turbulent wakes (see Results) were generated 
mor each test section configuration. The wake velocity 
defect and turbulence resulted from flow separation of the 


turbulent boundary layers on the cylinders. 


DEST SECTION 

As shown in Fig. 3, the test section consisted of two 
Eo llel, vertical side walls separated by flexible and 
EwUStable top and bottom walls, all made of Plexiglas. The 
flexible walls were made from one-eighth inch Plexiglas and 
fastened to the nozzle with piano hinges. Eight adjustmen-^ 
ponts, each with two threaded rods, allowed these walls tc 
be positioned to any desired area distribution and horizon 
dally levelled so that the test section geometry woulda re- 
main rectangular throughout. Appendix A contains the devel- 
opment of the relationship between the pressure gradient 
ENG area distribution, based on a one-dimensional analysis. 
Table I contains the dimensions necessary to ee area 
profiles and form the two pressure distributions used in 
es study. 

The right hand vertical wall of the test section contained 


thirteen slots through which the pitot-static probes were 
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inserted into the flow. These siots were nineteen inches 
high, three-eighths inch wide and spaced five inches apart. 
A teflon sliding seal prevented the influx of air through 
the slot in which the probe was being traversed. A filler 
plug for each unused slot was carefully hand milled and 
polished to make it conform perfectly with the flow side of 
the wall, thus eliminating the possibility of the plugs 
Protruding into the flow and inducing Step ahua t Yo rm Aithin 
aluminum strip wedged over each filler plug insured a tight 
B. 

The left hand vertical wall of the test η 
three rectangular access holes equally spaced along the 
axial centerline. These enabled the operator to remove and 
insert the wake generators, measure the probe position 
accurately, and level the lower adjustable wall. Also, along 
the axial centerline of this wall were thirty-three static 
pressure ports spaced one-and-a-half inches apart. A 
manometer bank connected to these ports gave a visual indi 
cation of the linearity and uniformity of the pressure dis- 
tribution. The manometer also served to indicate any regions 
of separation, as did several rows of tufts which were 
attached to the walls of the adverse pressure gradient test 


section. 


E. DIFFUSER AND MOTOR 
The test section was followed by two seven-inch square, 


four feet long steel sections, a diffuser section and a 
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centrifugal blower driven by an electric motor shown in 
Figure 4. The diffuser was attached to the blower by a 
rubber buffer section which reduced the amount of bivration 
transmitted from the motor to the tunnel. The centrifugal 
blower, operated at three thousand revolutions per minute, 
pulled air through the tunnel and exhausted it vertically 
through a muffler. The total mass flow through the tunnel 
was controlled by a damper mechanism just prior to the 


muffler. 


E. DATA SENSORS 

All data were taken with two calibrated Prandtl-type 
pitot-static probes, one measuring three-and-a-half inches 
mee che other eighty-three hundredths inch, stem to tip, 
thus allowing acquisition of data between test section data 
stations. Two U-tube manometers were connected to the probes 
ee hat direct readout of local static and dynamic pressure 
was obtained. These manometers, shown in Fig. 5, could be 
Melined at any angle between two and forty-five degrees, 
providing a wide range of sensitivity. A similar manometer 
was Ea to the nozzle static pressure ports so that the 
test section inlet dynamic pressure could be continuously 
Monitored. 

The traversing mechanism used to position the pitot- 
static probes in the wake was hand-operated and completely 
detachable from the tunnel. The probe was clamped firmly in 


its movable support fixture on the traverse which was moved 
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vertically by a screw drive mechanism. Probe position could 
be determined to within a hundredth of an inch with the aid 


of a ruled indicator on the traverse. 
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III. EXPERIMENTAL PROCEDURE 


A. TEST SECTION PREPARATION 

The test section was configured to produce both a 
favorable-only pressure gradient (Figure 6), and an adverse- 
favorable pressure gradient (Figure 3) over its entire 
length. For ease in setting the flexible walls to the de- 
Sired profile, two reference lines were scribed at the top 
and bottom of the right hand vertical wall of the test sec- 
tion. All dimensions given in Table I are referenced to 
Ese lines and indicate the distance in inches from these 
reference lines to the flow side of the flexible walls. 

The following procedure was used in setting the flexible 
EENNIS, Proceeding from the test section inlet, the right 
side of each adjustment fixture was set to the desired posi- 
on and tightened. After all eight fixtures had been sets 
each was re-checked to insure that its setting had not been 
Geecurbed. Then, "Т the aid of a bubble level, the left 
side of each adjustment fixture was moved until the wall was 
level. — the walls could onlysbewset precisel, 
to the desired dimensions at the eight adjustment points, 
the resulting area distribution in the test section never 
varied more than five percent from the desired area. The 
only variance occured when the axial area gradient was large, 
such as at the beginning of the favorable pressure gradient 


profile or in the transition section of the adverse-favorable 
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pressure gradient profile; otherwise the wall profile 


conformed exactly to the desired area distribution. 


B. DATA ACQUISITION 

All data were recorded from the water-filled U-tube 
manometers (Figure 5) connected to the pitot-static tubes 
and nozzles. A scale was placed between the tubes of each 
manometer providing a minimum reading of a hundredth of an 
inch deflection on the slant height. To eliminate any 
possible parallax and to assure consistent readings from 
time manometers, a portable plastic sight unit was placed 
mer the manometer scale. This unit restricted one's line 
of sight to a plane perpendicular to the manometer. 

At each position of the probe in the wake, the manometer 
readings were recorded and plotted on a graph so that it was 
immediately obvious when more data was needed in a particular 
portion of the wake. This technique was particularly advan- 
tageous in determining the exact center of the wake and the 
point where the freestream was encountered. 

The high turbulence intensity in the central core of the 
wake caused large fluctuations in the manometer water columns. 
These fluctuations were often as great as two-tenths inch 
even at manometer inclinations of thirty degrees. By ob- 
serving the fluctuations for several minutes at each probe 
position in the wake, it was possible to "average" these 
mhuctuations by eye and obtain a reading that was accurate 


to within an estimated five-hundredths of an inch. 
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The inclination (set with an inclinometer) of each of 
the three U-tube manometers was varied according to tunnel 
eBperating conditions and axial location of the pitot-static 
peobe. For a given axial data station and local flow veloc- 
ity, an inclination was chosen that would yield approximately 
Feten to eleven inch difference in» water eolumn slant heiche 
when the probe was positioned in the freestream. This 
assured adequate sensitivity and readability of the water 
columns at the lower velocities of the wake centerline. 

Traversing of the wake by the pitot-static probe was 
always done from the upper free stream edge to the lower 
free stream edge. The free stream was located by moving 
the probe away from the upper wall until the maximum slant 
height difference on the dynamic pressure manometer remained 
unchanged for a few tenths of an inch below the point at 
which it first occurred. This point was taken as the free 


Stream. The lower free stream was found in the same way. 
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IV. DISCUSSION OF RESULTS 


B TEST SECTION 

A Qualitative indication of the linearity of the tes: 
section pressure gradient in both the favorable and the 
adverse-favorable configurations can be seen in Figures 7, 
8. The limitations of the test section in generating an 
Exactly linear pressure distribution result from the lack 
of completely flexible upper and lower test section wall 
and from the limited number of wall adjustment points from 
which to set the area profile. Considering these physical 
imitations and the fact that the desired area distribution 
in the test section was derived from one-dimensional analysis, 
Bue resulting actual pressure distribution in the test 
Section is quite good. A comprehensive program to calibrate 
the test section was not undertaken since only a qualitative 


idea of its performance was necessary for this investigation. 


Eee EXPERIMENT 

The experiment was designed to compare the decay of two 
Similar wakes which have had different flow histories. This 
was done, within the physical constraints of the test section, 
by using a two-inch diameter wake generator in a favorable 
pressure gradient channel sixty eight inches long and a one- 
inch diameter wake generator in an adverse-favorable pressure 


gradient channel with the favorable pressure gradient 
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occupying the downstream thirty four inches (Figures 3 and 
8). The two wakes are compared within the favorable 
pressure gradient regions downstream of the point where 

both wakes have the same velocity defect. The approximate 
Aution of the integral momentum equation in Appendix C 
illustrates the possibility that an adverse pressure gradient 
can sustain a wake deficit, such as the one studied here. 

fm addition, the wakes are compared at the same relative non- 
dimensional distance into the favorable pressure gradient. 
This was possible because the favorable pressure gradient 
Menpths were in the same ratio as the initial wake widths. 
It is apparent from the experimental data in Figure 9 that 
the velccity defect decays rapidly when the flow enters the 
B Uorable pressure gradient directly from the nozzle. THONS 
ever, the velocity defect persists [Figure 10] and then 
decays at a much slower rate [Figure 11] when the flow has 
Previously encountered an adverse pressure gradient. Since. 
mee rate at which the velocity defect diminishes is direct у 
related to the eddy viscosity, it can be at least qualita- 
tively concluded from Figures 12,13 that the eddy viscosity 
EE turbulent wake is a function of the flow history and 
mes effect on the turbulent field. Classically, it is 
@esumed that the eddy viscosity is purely a function of the 
local mean flow properties, we must conclude that care must 
be exercised when applying classical eddy viscosity models 


to pressure gradient situations. In order to test this 
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conclusion, a numerical analysis of the experimental data 


Was conducted. 


C. COMPUTER PROGRAM 

The data analysis computer program uses a forward 
difference scheme to solve the equations of motion for the 
turbulent shear coefficients (CTAU) and the eddy viscosi- 
ties (EDDY) as described in Appendices B and D. The exper- 
imental data obtained from the wind tunnel were used to 
test the program and assure correct processing with no 
Eng errors. The numbers generated for CTAU and EDDY were 
unreasonably large and therefore the results are inconclusive. 
fers felt that the increments in the data in the axial 
direction used in finite differencing the data were too 
Targe for an accurate approximation of the partial deriva- 
tives which appear in the equations (see Appendix B). A 
better method would be to smooth the data and use a central 
difference numerical scheme to approximate the partial 
derivatives. The data taken in the experimental ΗΕ σε 
this investigation were not adaptable to the central 
difference bonne. Given experimental data with sufficiently 
small increments in the axial direction, the program will 
Provide the correct values for the parameters listed in 


Appendix D. 
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V. CONCLUSIONS 


The experiment described in this report was oec mered со 
test the validity of classical eddy viscosity models in 
free turbulent mixing situations with axial pressure gra- 
ПЕ nts, The results indicate that those elassrealsmegels 
Dich are explicitly independent of the local turbu@ence 
structure will fail to predict the experimentally observed 
pmeemomena which are a strong function of the turbulent field. 
In these pressure gradient situations, the local turbulent 
kinetic energy is not a unique function of the local non- 
Aw bulent mean flow and the influence of the turbulence 
on the form of the Reynolds stresses must be an explicit 


Ome En any eddy viscosity model. 
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Figure 1 


View of tunnel showing nozzles used with (1) fav- 
orable pressure gradient and (2) adverse-favorable 


pressúre gradient test section σας Еа оа: е 
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Figure 2 


View of Plexiglass nozzle and test section showing 
(1) 2-inch wake generator installed in the nozzle 


and (2) the pitot tube traverse mechanism. 
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Figure 3 


View of Plexiglass test 
adjustable walls in the 
gradient configuration, 
these walls and (3) the 


tubes are inserted. 


section showing (1) the 
adverse-favorable pressure 
(2) the adjusting rods for 


slots through which the pitot 
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Figure 4 


Overall view of the wind tunnel showing the diffuser, 


centrifugal blower and its muffler at extreme right 


"ου. Πεντε. 
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Figure 5 


View of inclined U-tube manometers used to obtain 
pressure readings in the nozzle and test section. 
The sight unit (1) used to eliminate parallax in 


reading the manometers is shown on the left. 
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Figure 6 


Test section configured for the favorable pressure 
gradient. The milled slot through απο του 


tube is inserted can be seen at station H. 
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Figure 7 


View of manometer bank connected to the static 
pressure ports along the test section centerline. 
Shown is the linear pressure distribution in the 
test section resulting from the configuration in 
Fig. 6. The horizontal water columns are not 


connected to the tunnel. 
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Figure 8 


View of manometer bank connected to the static 
pressure ports along the test section centerline. 
Shown is the pressure distribution in the test 
section resulting fem the соп тепаа С 7 
The horizontal water columns are not connected to 


the tunnel. 
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APPENDIX A 


DESIGN OF THE TEST SECTION AREA DISTRIBUITON 


The local cross-sectional area of спе ΠΠ. 
the static pressure of the fluid flowing through the section 
are related through the equation of continuity and Bernoulli's 
equation. For the steady, isentropic flow of an incompres- 


prole fluid, these equations can be written: 
ПА = constant Gir) 
р + ος = constant | (2) 


2 


In the log-derivative plane (1) and (2) become 


av _ алд 
se x (3) 
dp = a = (4) 


Substitution of (3) into (4) yields 


ма (5) 


dp = pV А 


From (1) V = C4/A, where Ci is the volumetric flow rate at 


the inlet to the test section. Thus, (5) can be written 
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After specifying a particular pressure gradient, (6) can 

be solved to determine the required area distribution. 
For a constant pressure gradient, ср = -C,, and inte- 

gration of (6) for this simple case is straightforward. 


The result can be written 


Alx) il 
MS AA (7) 


T 1 + Bx 





where ß = C,/9, >» А, is the test section inlet area and q; 
is the dynamic pressure at the test section inlet. For a 


Specified inlet and exit area 
| 2 


umere L is the length of the pressure gradient, and A, is 


the channel area at x = L. 


The area profiles for the favorable and adverse-favorable 


pressure gradients used during this investigation are shown 


in Figs. 3 and 6, respectively. Table I contains the specific 


Gimensions used to set the flexible walls of the test section 


to the desired area profile. Shown for comparison in this 


table are measurements of the upper and lower flexible walls 


as they were set prior to taking data. The differences in 
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dimensions, where they exist, between the theoretical and 
actual wall settings show the degree of conformity of the 


flexible walls to the desired profile. 
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APPENDIX B 
DEVELOPMENT OF AN EQUATION FOR CALCULATION 
OF THE LOCAL TURBULENT SHEAR COEFFICIENT 
The momentum equation governing the mean flow properties 
of a steady, two-dimensional, incompressible turbulent 
shear flow can be written, using Reynolds' averaging 


technique [Ref. 3]: 


— gu = gu 1 dp l дт 
== = Z — — >n t 268 
М 9х БЕ ду p dx PO y = 
where t = -pu'v' and suis has been neglected with respect 


to p. The same arguments which Prandtl used to derive 
equation (1) from the Nee equations, in the case 
of a boundary layer, are directly applicable to a two- 
dimensional wake in approximately parallel flow through a 
ENCLangular duct.. Introducing the integrated form of the 


continuity equation 


y — 
AES К ES dy (2) 


Ento equation (1) transforms it to 


— ο y — — 
qu u ou = ὍΣ 


У о 


By non-dimensionalizing equation (3) with respect to the 


local freestream velocity, Us» and utilizing the τοι; 
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of differential calculus, equation (3) becomes 





u u 
l u πιο x u Y ú u С = 
— (—), + (—)” — - (—) П (=) dg = 288 (u 
De Eo 2 Ho Mo Че y 0 us 5 ча У Ue 0 u 
@ 
2 1 P + n (4) 
PU, 


where C. = 21/pu ° and subscripts x and y indicate partial 
differentiation with respect to x and y. Equation (4) is 


now integrated from y = 0 toy = y an arbitrary point in 


ct 
the wake, and relax the condition that equation (!) be 
satisfied at every point in the flow, allowing it to be 
satisfied on the average between y = 0 and y = y After 
the necessary integration by parts and considerable manipu- 
lation, the final equation for the turbulent shear coeffi- 


cient, ar can be written 











guy ) In = u(y_) — e Ут: 
ET u m u u x l- = 
= f [2(—3 - JE). + (Jay + / Р ау 
2 0 Un Un u, X U, u. 0 pu, X 


(5) 


Equation (5) was solved numerically with the computer 
program listed in Appendix D. 


The validity of equation (5) can be checked by letting 


ME 
Ui ye) eu δαν же 
е? ον m 
aP a - 
x E e dx 





Thus PEE) - 0 and equation (5) reduces to 


du 
de AS x` = 
àx*u ду 20 ThE 0 (6) 


where 


Equation (6) is the well known momentum integral equation 


in boundary layer theory. 
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APPENDIX C 
SOLUTION TO THE MOMENTUM 
INTEGRAL EQUATION 
The momentum integral equation obtained in Appendix B 
can be integrated in the special case of a constant pressure 


gradient. From Eq. (6) in Appendix B 


du 


-- -ᾱ (294 6*) = 0 (1) 


u 
e 


+ 


we 
> | 


where 6, 6* are the momentum and displacement thickness 


defined in Appendix B. 


du 
An expression for dr сап be obtained from the inviscid 
Bernoulli equation 
ve. em (2 
dx pu ах 


Substitution of Eq. (2) into Eq. (1) yields a differential 


equation in terms of 9, m and H = 8%*/0 (the shape factor). 








ax 
q 1 dp w 
a, (In 0) πα ο ο. (3) 
pu 
e 
The solution of Eq. (3) is 
0 Ὁ μι KC: 
= = ех / (SP) (H + 2)dx] (4) 
8 2 `dx 
1 0 pu, 
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From Appendix A, T - -Ва+ and the integrand of Eq. (4) 


can be written 


s. (MA) (b) 
оп 5 


From the results of Appendix A, the ratio is a known function 


or Xx 


© 
= 


_ 1 
a. TITR n 


O 
Ф 


After substituting Eqs. (5), (6) into Eq. (4) and integrating, 


the final solution becomes 


_ 1 а BH 
$7 TGUER a 5{т + ху αχ] (1) 


D 


Thus, experimentally derived values for the shape factor 

@an be substituted into Eq. (7), and the resulting value of 

8 can be compared with the experimentally obtained values. 
For the restrictive case, H = constant, we can extend 


Beis solution to 


Ө 1 
RU ce (8) 
θη (3 « gx) R72 


Further, for H = 1.0 (a close approximation under the 


conditions of this investigation), Eq. (7) reduces to 


HO 





ug 1 (9) 
θη (1 + вх) 3/2 


This result compares favorable with experimental data апа 
serves as a guide for calculating the growth or decay of 


wakes in pressure gradients. 
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APPENDIX D 


DATA REDUCTION PROGRAM 


A. GENERAL 

The following FORTRAN IV computer program uses the 
experimental values of dynamic and static pressure through 
the wake in inches of water and calculates the wake momentum 
thickness, displacement thickness, turbulent shear coeffi- 
cient, eddy viscosity, and Prandtl mixing length. ~The 
results are printed out as well as plotted on the calcomp 


Potter. 


B. NOMENCLATURE 

The following parameters and arrays are used in the 
program, with their respective dimensions shown in 
parentheses: 


ALV - array used to store interpolated values of 
dynamic pressures (feet) 


AMIXL - local Prandtl mixing length in the wake (feet) 


ΕΠΙ = local turbulent shear coefficient for each 
point in the wake 


DELRAT - ratio of succeeding DELSTR's to that at the 
first data station 


DELSTR - wake displacement thickness (feet) 


DELX - axial distance between succeeding probe tip 
positions (feet) 


DIAMTR - diameter of wake generator (inches) 


EDDY - eddy viscosity 


D 





ВЕЗЕТ, 
LABEL 
MN 


MX 
MY 


PATM 
Q 
SHAPE 
THETA 
THETAP 


THETAV 
THRAT 


ЕЕЕ 
TITLE4 


U 
URATIO 
XIN 
XOUT 


XPOS 


residual in momentum integral equation when 
experimental values for THETA and SHAPE are 
Sulbar L vtec Waro Lt 

label for subroutine DRAW 


number of data points laterally through the 
wake 


number of axial data positions 


number of points (always odd) into which the 
wake is divided 


array of local static pressures for each point 
in the wake (psf) 


atmospheric pressure (inches of mercury) 
freestream dynamic pressure (psf) 

wake shape factor (DELSTR/THETA) 

wake momentum thickness (feet) 

inclination of static pressure manometer (degrees) 


inclination of dynamic pressure manometer 
(degrees) 


ratio of succeeding THETA !s to that at the 
first data station 


title tor subroutine DRAW 
title for subroutine DRAW 


array of velocitiesifor each polo merke 
wake (feet per second) 


ratio of local wake velocity to local free- 
stream velocity 


array of probe traverse positions through the 
wake (inches) 


interpolated probe traverse positions for each 
of the MY points in the wake (inches) 


probe tip position measured from downstream 
side of the wake generator (inches) 
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YOUT - interpolated manometer readings (for YQ and 
YPS) for each of the MY points through the 
wake (inches) 


τες - static pressure manometer slant distance 
height difference (inches) 


YQ - dynamic pressure manometer slant distance 
height difference (inches) 


SUBROUTINES 


INTEG - uses Simpson's rule to numerically integrate 
the function F. It must be опере еи 
points (always odd) equally spaced, DELX 
units. apart. 


PIF2 ~ second-order interpolation function. X is the 
point for which interpolation is desired. 
XLIST and FLIST contain the N number of abscissas 
and ordinates, respectively. one ciemreo le re 
be interpolated. 


DECIDE - insures that the wake is divided into an odd 
number of points spaced 0.1 3meBes apart. ΙΤ 
the wake contains an even number of experimental 
data points (found by subtracting the lower 
traverse position from the upper position and 
multiplying by ten), ОЕСР K ο-ᾱ ου n 
the upper edge of the wake and sets its value 
equal to the freestream value. The total 
number of points are printed out. 


FILL - insures that the wake at succeeding downstream 
data positions contains as many data points as 
the one before it. πας ae cols 
precluded in the computation scheme used for 
EPSIL and CTAU. If the succeeding data station 
has fewer points in its wake, FILL copies the 
uppermost wake point (the freestream point) 
into additional points as necessary to equal the 
number of points@wat Te zeUrzentspesition. 


DRAW - SSP3 subroutine for the IBM 360/67 at the Naval 


Postgraduate School which uses the CALCOMP 
plotter to graph thezpreenenszeutpür. 
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D. USAGE 

The program will accept up to twelve data positions in 
the streamwise direction. If more than that number are 
required, the dimensions of the appropriate arrays (ALP, 
pve U,P,A,B,C,D,EDDY, AMIXL, XOUT, URATIO, CTAU Εμ ТИЕ Б 
SHAPE, EPSIL,Q,DELRAT, THRAT, THETAV , THETAP, PAD ees ΡΕ. 
NUM) will have to be modified accordingly. The program 
is set up for exactly three complete runs (a run is defined 
as the processing of data taken for a particular test 
section configuration and inlet dynamic pressure). If more 
or fewer runs are necessary, the limits for the KK DO-LOOP 
have to be modified accordingly. There is no restriction 
on the number of data points that can be obtained at each 
axial data position. Subscript I always refers to each 
axial position for which data was taken, and subscript J 
always refers to the points through the wake, at a given 
axial data position. 

The following E data is supplied by the user in the 


order presented. Input format is given in parentheses. 


wo cards: TITLE] (6A8) - first graph title 
Two cards: TITLE! (6A8) - second graph title 


One card: LABEL (12A4) - curve label. One number in each 
four columns to designate the data position, 
i.c. 1-2 3 "ELO. 


For each run the following data is needed: 


One card: run identification (up to спасе сеге) 
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One 


опе 
Опе 
Опе 
Опе 
Опе 
Опе 
Опе 


Опе 


сага: 


ог тоге 
ог тоге 
or more 
or more 
card: 

or more 
or more 


or more 


PROGRAM 


number of data positions (I2); wake generator 
diameter, DIAMTR 'inehes , (F10.2) 


cards:  THETAV, degrees (8F10.1) 

cards:  THETAP, degrees (8F10.1) 

cards: PATM, inches (8F10.2) 

cards: XPOS, inches (8F10.2) 

data position identification (up to 47 ers) 
cards: XIN, inches (8F10.3) 

cards: YQ, inches (8F10. 3) 


cares: YPS, inches (8F10.3) 


LISTING 


A listing of the program which has been described here 


is included on the following six pages. 
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